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ABSTRACT 

We present a detailed analysis of the number count and photometric redshift distribution of faint 
galaxies in the Hubble Deep Field (HDF), paying a special attention to the selection effects including 
' the cosmological dimming of surface brightness of galaxies, under the observational condition employed 

. in this field. We find a considerably different result from previous studies ignoring the selection effects, 

(N ■ 

and these effects should therefore be taken into account in the analysis. We find that the model of 
'"^ ■ pure luminosity evolution (PLE) of galaxies in the Einstcin-de Sitter (EdS) universe predicts much 

\ smaller counts than those observed at faint magnitude limits by a factor of more than 10, so that a 

' very strong number evolution of galaxies with rj > 3-4 must be invoked to reproduce the /814 counts, 

when parametrized as cx (1 + z)''. However we show that such a strong number evolution under 
■ realistic merging processes of galaxies can not explain the steep slope of the i?45o and Vgoe counts, and 

^ . it is seriously inconsistent with their photometric redshift distribution. We find that these difficulties 

\^ I still persist in an open universe with f2o ^ 0.2, but are resolved only when we invoke a A-dominated 

CN| ' flat universe, after examining various systematic uncertainties in modeling the formation and evolution 

JnI I of galaxies. The present analysis revitalizes the practice of using faint number counts as an important 

' cosmological test, giving one of the arguments against the EdS universe and suggests acceleration of the 

, cosmic expansion by vacuum energy density. While a modest number evolution of galaxies with 77 < I is 

' still necessary even in a A-dominated universe, a stronger number evolution with 77 > I is rejected from 

the HDF data, giving a strong constraint on the merger history of galaxies. 

^ ■ Subject headings: cosmology: observations — galaxies: evolution — galaxies: formation 

^ ; 1. INTRODUCTION Takahara 1988; Fukugita, Takahara, Yamashita & Yoshu 

• • AT u 4-- r c ■ 4- T ■ ■ c 4-1, 4- 1990: Yoshii & Peterson 1991). However, the obtained 

J> Number countmg of tamt galaxies is one of the most . . ' ' 

.—I r J i 1 u 4.- ^ J. J. ■4-^ i • i ii r Constraints On cosmological parameters bascd On the PLE 

j/Tj ' fundamental observational tests with which the forma- ° ^ 

: tion/evolution of galaxies as well as the geometry of the "^"^^^^ ^^^^ '^^^ ^^^^ deterministic because of pos- 

3: . . ujrpii,- 4. J 4- c 4.U sible number evolution of galaxies by mergers. Particu- 

universe is probed. Ihe best view to date of the optical ° j & 

r-.n 1 ,• . uj-i,tj uui T7- ij larW, whcu thc PLE model is used, the Einstein-de Sit- 

sky to tamt nux levels is given by the Hubble Deep i'leld ' 

(HDF, WiUiams et al. 1996), and it provides a valuable ^^'^^^ ^^^^""'^ = underpredicts the observed 

information to a wide range of studies on galaxies and S^^^^^ ^^'""^ magnitudes, but a simple model 

, A u • J- J f 4-u tjTM? 1 of galaxy number evolution can reproduce the observed 

cosmology. A comprehensive study of the HDl galaxy t> j t- 

,1, r jun ii-ii nnno^ j counts and save the EdS universe (Rocca-Volmerange & 

counts has been performed by Pozzetti et al. (1998), and ^ ° 

4-u f A i-u 4- ■ ^ A ^ f \ ■ -4- 1 Guiderdoni 1990; Pozzetti et al. 1996). This degeneracy 

they found that a simple model or pure luminosity evolu- ' j o j 

,. /-DT t:^\ ■ i,-v,i- 1 -ijj-i between the effects of galaxy evolution and cosmology has 

tion (PLE), m which galaxies evolve passively due to star . 

r ^ ■ , ■ ■,^ , 1 1 , • been a maior problem when one uses the galaxy number 

formation histories without mergers or number evolution, j t- o j 

gives a reasonable fit to the HDF counts in all the four ^^^^^^ ^° determine the geometry of the universe. 

\- J c TT T, Tr J T ^ • -ii The information of redshifts is able to break such a de- 

pass bands of (7, £>, V , and i, when an open universe with 

i7 — 1 is assumed generacy, because luminous galaxies at great distance are 

mi ■ r , 1 1 r 1 ■ -ii ii ■ distinguishable from dwarf galaxies in a local universe. Al- 

Ihe increase of the number of galaxies with their ap- ° ° 

, ., J • • 11 J c though most of the HDF galaxies are too faint to mea- 

parent magnitude was originally proposed as a measure of " ° 

,1 , r , , • /c J 1nc1^ J -J sure the spectroscopic redshifts, several catalogs of their 

the geometry of the universe (Sandage 1961), and consid- t- t- j e> 

erable efforts have been made along this line (e.g., Yoshii & 
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photometric redshifts have been pubUshed (Sawicki, Lin, 
& Yee 1997; Wang, BahcaU, & Turner 1998; Fernandez- 
Soto, Lanzetta, & Yahil 1999). The follow-up studies 
based on these catalogs show that the photometric red- 
shifts give reasonably reliable estimates of spectroscopic 
redshifts and are useful for a statistical study of the HDF 
galaxies. Here we give a combined analysis for the HDF 
counts and redshifts and constrain the cosniological pa- 
rameters separately from the merger history of galaxies. 

Both the number count of faint galaxies and their red- 
shift distribution are significantly affected by the selection 
effects inherent in the method of detecting galaxies in faint 
surveys, but these important effects have been ignored in 
almost all previous studies except for Yoshii & Fukugita 
(1991) and Yoshii (1993). It is well known that the sur- 
face brightness of galaxies rapidly becomes dimmer with 
increasing redshift as c>c (l + z)-'^ (Tolman 1934), and this 
cosmological dimming makes many high-redshift galaxies 
remain undetected below the threshold value of surface 
brightness adopted in a galaxy survey (Pritchet & Kline 
1981; Tyson 1984; EUis, Sievers, & Perry 1984). The see- 
ing or smoothing of an image furthermore lowers its sur- 
face brightness, and the photometry scheme used in a sur- 
vey heavily affects a magnitude estimate of the faintest 
galaxies. Some observers apply corrections to raw counts 
of faint galaxies for those undetected, but it is in princi- 
ple difficult and heavily model-dependent to estimate the 
number of undetected galaxies. Rather, the best way is 
to make theoretical predictions with the selection effects 
taken into account and then compare them directly to raw 
counts (Yoshii 1993). 

This paper is the first analysis of the HDF galaxies in 
which the above selection bias against high-redshift galax- 
ies is explicitly incorporated. We use a standard PLE 
model of galaxies including the effects of internal dust ob- 
scuration and intergalactic HI absorption. Number evo- 
lution of galaxies is also allowed for with simple modifi- 
cations to the PLE model. Throughout this paper, we 
use the AB photometry system with the notation of U300, 
B4501 Veoe, and /814 (Williams et al. 1996). In we 
present a detailed description for models of galaxy evo- 
lution and formulations to calculate galaxy counts and 
redshift distribution with the selection effects taken into 
account. Extensive calculations of number count predic- 
tions and comparison to the HDF counts are given in §|[ 
checking in great detail the uncertainties arising from the 
prescribed properties of local galaxies and their evolution. 
We will give the comparison of the model predictions with 
the observed photometric redshift distribution in We 
discuss the results in §|^. The summary and conclusion of 
this paper are given in §^ 

2. THE MODEL OF GALAXIES AND DETECTION 
IN THE HDF 



First we describe the basic ingredients involved in our 
theoretical modeling such as the local luminosity function, 
galaxy evolution in luminosity and number, internal and 
intergalactic absorption, and the selection effects. Then 
we will present the formulations to calculate the number 
count of faint galaxies and their redshift distribution. 

2.1. Galaxies at Present, and Their Evolution 

We use a standard PLE model of galaxy evolution, in 
which galaxies are classified into five morphological types 
of E/SO, Sab, Sbc, Scd, and Sdm. Spectral energy distri- 
butions (SEDs) and their evolution are calculated by using 
the galaxy evolution model of Arimoto & Yoshii (1987) 
for elliptical galaxies and the II model of Arimoto, Yoshii, 
& Takahara (1992) for spiral galaxies. These models are 
constructed to reproduce the photometric and chemical 
properties of present-day galaxies. In order to see the 
systematic uncertainty in evolution models, we will also 
use an updated version of these models by Kobayashi et 
al. (1999) using the latest database of stellar populations 
compiled by Kodama & Arimoto (1997). We set the epoch 
of galaxy formation at z^? = 5 as a standard and change 
this value to see the systematic uncertainty. 

The luminosity function of local galaxies is also impor- 
tant in predicting the number count and redshift distribu- 
tion of faint galaxies. We use the type-dependent (E/SO, 
Spiral, and Irr) _B-band luminosity function derived from 
the Second Southern Sky Redshift Survey (SSRS2, Marzke 
et al. 1998). We associate the Sab, Sbc and Scd models 
to be assigned to spiral galaxies, whereas the Sdm model 
assigned to irregular galaxies. The relative proportions of 
Sab, Sbc, and Scd are taken from Pence (1976). In order 
to check the systematic uncertainty related to the luminos- 
ity function, we also use the type-independent luminosity 
function of Stromlo-APM redshift survey (Loveday et al. 
1992) and the type-dependent luminosity function from 
the Center for Astrophysics (CfA) redshift survey (Huchra 
et al. 1983). Their Schechter parameters are tabulated in 
Table |l] (see also Efstathiou, Ellis, & Peterson 1988 and 
Yoshii & Takahara 1989). 

2.2. Absorption 

The above models of galaxy evolution do not include 
the absorption by interstellar dust which becomes signif- 
icant for high-z galaxies when observed in optical bands. 
In order to take this effect into account, we make a phys- 
ically natural assumption that the dust optical depth is 
proportional to the column density and metallicity of the 
gas. In fact, it is well known that the Galactic extinction 
is well correlated to the column density of the HI gas (e.g., 
Burstein & Heiles 1982). It is also known that the dust 
opacity becomes smaller in order of decreasing metallic- 
ity from the Galaxy to the Large and then Small Magel- 
lanic Clouds, when the gas column density is fixed (e.g., 
Pei 1992). Since the galaxy evolution models give the gas 
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fraction fg and the metallicity Zg in the gas, the dust opti- 
cal depth is calculated from r^ust = i'^fgZgr~'^{M / Lb)Lb, 
where rg, M, and Lq are the effective radius, the baryon 
mass, and the _B-band luminosity of a galaxy, respectively. 



(We will describe the treatment of galaxy size in §2.3.) The 
proportionality constant k is chosen to be consistent with 
the present-day, average extinction of Ay ^ 0.17 taken 
from the Galactic extinction map (Burstein & Heiles 1982; 
Schlegel, Finkbeiner, & Davis 1998) and a theoretical es- 
timate (Hatano, Branch, & Deaton 1998). The standard 
extinction curve of our Galaxy (e.g., Pei 1992) is used for 
the wavelength dependence of the optical depth. 

Given the optical depth, the attenuation of emerging 
stellar lights depends on the spatial dust distribution. Fol- 
lowing Disney, Davies, & Phillipps (1989), there are two 
extreme cases such as the screen model in which the dust 
is distributed on the line of sight to stars, and the slab 
model in which the dust has the same distribution with 
stars. Neglecting the scattering of lights by dust, the at- 
tenuation factor of stellar lights is given by exp(— Tdust) for 
the screen model and [1 — exp(— rcjust)]/''dust for the slab 
model. In fact, the galaxy evolution model used here has 
been made to reproduce the present-day SED of galaxies 
which has been already affected by dust obscuration. We 
take into account this point and hence correct the above 
attenuation factors by using the optical depth at present. 
Therefore the present-day SEDs of model galaxies are the 
same for all the prescriptions of dust-free, screen, and slab 
models. 

In the slab model the apparent reddening reaches an 
asymptote when the optical depth becomes much larger 
than unity, because the observed lights are emitted from 
surface regions of a galaxy where the optical depth to an 
observer is low. However, the observed correlation be- 
tween the power-law index of UV spectra and the Balmer 
line ratio, both of which are a reddening indicator, ex- 
tends well beyond the asymptote. This indicates that the 
observed reddening of starburst galaxies is larger than ex- 
pected from the slab model, and at least some fraction of 
dust should behave like a screen (for detail see Calzetti, 
Kinney, & Storchi-Bergmann 1994). 

We then use the screen model as a standard, consider- 
ing that UV and optical observations of starburst galaxies 
favor the screen dust. It may be an extreme prescription 
that all dust is distributed as a screen, but note that the 
emergent lights are rapidly attenuated exponentially once 
a considerable fraction of dust contributes to the screen. 
Therefore the screen model is more appropriate than the 
slab model in which the attenuation factor decreases only 
moderately like T^g^ when Tdust — > 

The spectral energy distributions (SEDs) of galaxies of 
various types are given in Fig. ^ at several epochs of galaxy 
evolution for three cases of dust extinction such as the 
screen model (solid line), the slab model (dashed line), 
and the no-extinction model (dotted line). The effect of 



dust extinction is especially important for elliptical galax- 
ies at high redshifts, where UV radiation is quite strong 
because of the initial starbursts supposed in the galactic 
wind model of elliptical galaxies (Arimoto & Yoshii 1987). 
According to the method described above, the dust optical 
depth during the initial starburst phase of elliptical galax- 
ies is estimated to be much larger than unity (r > 10) for 
UV photons at ~ 2000 A in the restframe. On the other 
hand, the evolution of dust obscuration makes UV lumi- 
nosity of late-type spiral galaxies brighter at early epochs 
than that without the dust effect, because of the lower 
metal abundance than the present-day galaxies. However, 
this effect in late-type spiral galaxies is not significant be- 
cause they are not heavily obscured at present. 

In addition to the internal absorption by dust, the in- 
tergalactic absorption significantly affects the visibility of 
high-rcdshift galaxies (Yoshii & Peterson 1994; Madau 
1995). Lights from a distant galaxy at rest wavelengths 
below the Lyman limit (912 A) and those below the Ly- 
man a line (1216 A) are extinguished by Lyman continuum 
absorption and Lyman series line absorption, respectively, 
in intergalactic HI clouds along the line of sight. We in- 
clude this effect consistently in our theoretical calculations 
making use of the intergalactic optical depth calculated by 
Yoshii & Peterson (1994). The optical depth of this ab- 
sorption is shown in Fig. |2[ as a function of observed 
wavelength for various source redshifts. 

2.3. Selection Effects 

Apparent surface brightness and size of an image in a 
survey observation are the essential quantities for it to be 
detected as a real galaxy. We calculate these quantities 
of a model galaxy assuming its intrinsic luminosity profile 
and size and taking into account the cosmological dim- 
ming and the observational seeing. For the details of the 
formulations, see Yoshii (1993). 

2.3.1. Galaxy Sizes and Luminosity Profile 

In our PLE model, we assume that the galaxy size does 
not evolve except for the case of mergers of galaxies, and 
use the empirical relation between the effective radius 
and absolute luminosity Lb for local galaxies. If we allow 
for the number evolution of galaxies, we must take into ac- 
count the change of galaxy sizes, and we will discuss this 
in §2.4, Fig. ^ shows the size-luminosity relation of local 
elliptical and spiral galaxies. The data of elliptical galax- 
ies are taken from Bender et al. (1992), and those of spiral 
galaxies from Impey et al. (1996). Although there is a sig- 
nificant scatter in this relation, we use a simple power-law 
relation for this relation as 



oc L 



2.5/p 



(1) 



or, if expressed in terms of the absolute B magnitude, 
-Ms =plogre + g-H(p-5)log(iJo/50km/s/Mpc) . (2) 
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Elliptical galaxies form two distinct families which follow 
the well-separated sequences at low luminosities in the r^- 
Lb diagram. One is the ordinary sequence from giant 
through dwarf elliptical (GDE) galaxies, while the other is 
the bright sequence from giant through compact elliptical 
(GCE) galaxies (see Fig^). Since the predictions of galaxy 
number count are not sensitive to whichever sequence is 
used in the analysis (see §3.2 and Fig m, we take the GDE 



sequence as the standard size-luminosity relation of ellip- 
tical galaxies in this paper. 

The Ts-Lb relations fitted to the data in Fig.^j yield 
(p, q) = (6.0, 16) and (3.5,18.7) for the GDE and GCE se- 
quences of elliptical galaxies, respectively, and (9.4, 12) for 
spiral galaxies. In order to examine the uncertainty due to 
the significant scatter in the t^-Lb relation, we derive the 
standard deviation in A(logrg) from the best-fit relation 
and repeat calculations with the shifted relations shown in 
Fig. ^ by the dashed lines. 

The radial distribution of surface brightness is as- 
sumed to follow de Vaucouleurs' (1962) profile [S oc 
exp[— (r/re)"'^/'*] ) for elliptical galaxies and an exponen- 
tial profile for spiral galaxies (Freeman 1970). Then we 
can calculate the radial distribution of surface brightness 
of a galaxy at given redshift in any passband, when the 
galaxy type, the present-day B luminosity, and the evo- 
lution model are specified. This surface brightness profile 
should be convolved with a Gaussian point-spread function 
(PSF) having the same dispersion with the observational 
seeing. 

2.3.2. Object Detection 

Let ^th be the surface brightness threshold adopted in a 
galaxy survey. Strictly, this threshold could change due to 
different noise levels within a survey field, but we use a sin- 
gle value for the simplicity. When the observed profile of 
a galaxy image is calculated as above, we can estimate the 
isophotal size which encircles a region of a galaxy image 
with surface brightness brighter than Sth_- If the isophotal 
size of a galaxy is zero, i.e., the central surface bright- 
ness is fainter than 5'th, this galaxy can not be detected in 
the survey. Usually a minimum isophotal diameter Dmin, 
which is comparable to the seeing size, is adopted as a 
condition for an image to be detected as a galaxy. We can 
calculate the isophotal diameter for a model galaxy, then 
it is easy to check whether this galaxy meets the detection 
criterion in the galaxy survey. 

2.3.3. Photometry Scheme 

There are three photometry schemes to evaluate the ap- 
parent magnitude of galaxies such as isophotal, aperture, 
and pseudo-total magnitudes. These magnitudes may sig- 
nificantly differ especially near the detection limit, and this 
difference should be included in the theoretical modeling 
of the selection effects inherent in the method of detecting 
faint galaxies. 



The isophotal magnitude is the flux within the isopho- 
tal size of a galaxy image. The aperture magnitude is 
the flux within a fixed aperture which is adopted by ob- 
servers. Some observers often make corrections to these 
magnitudes into pseudo-total magnitudes, which are in- 
tended to mimic an ideal total flux of a galaxy. However, 
such corrections use a model luminosity proflle to evalu- 
ate the flux from an 'undetected' part of a galaxy image. 
It is inappropriate to make a model-dependent correction 
to the observed quantities with which various theoretical 
models are compared. The best way is, on the contrary, 
to incorporate all the necessary corrections in the theo- 
retical models to be compared directly with the observed 
quantities. We therefore suggest that observers should also 
present raw counts and magnitudes, in addition to present- 
ing corrected quantities in their papers. 



2.3.4. Detection and Photometry of the HDF 

Galaxies 

Here we describe the detection processes of the HDF 
galaxies, following Williams et al. (1996), and see also 
Bouwens, Broadhurst, & Silk (1998). Object detection is 
performed in the combined Veoe + ^8i4 image. It is flrst 
convolved with a flxed smoothing kernel of 25 pixels (=0.04 
square arcseconds), then pixels having values higher than 
a fixed threshold above a local sky background are marked 
as potentially being part of an object. After thresholding, 
regions consisting of more than contiguous 25 pixels are 
counted as sources. It corresponds to an isophotal diame- 
ter limit of Dmin ^ 0.2 arcsec, about 1.6 times larger than 
the FWHM of the PSF. 

The surface brightness threshold is not clearly indicated 
in Williams et al. (1996), but we can evaluate this value 
from the surface brightness distribution of the HDF galax- 
ies. In Fig. ^ we plot their apparent magnitudes ver- 
sus average surface brightness in the four passbands of 
C^300i ^450, ^606, and /814 by using the published sizes 
and magnitudes of the HDF catalog. Here the size refers 
to the isophotal size of the combined Veoe + ^814 image, 
and the magnitude refers to the isophotal magnitude mea- 
sured within the isophotal size. First we consider the Vgoe 
and /814 bands. Figure ^ shows that no galaxies are de- 
tected when the average surface brightness is fainter than 
St\i — 27.5 mag arcsec"^ in Veoe and 27.0 mag arcsec"^ in 
/814- The faintest surface brightness detected in the HDF 
for galaxies with a fixed isophotal magnitude becomes 
brighter for brighter galaxies. This trend occurs because 
the size of galaxies is larger for brighter galaxies. The 
edge of an isophotal image corresponds to the isophotal 
limit, and hence its average surface brightness within the 
isophotal area is always brighter than its threshold when 
the galaxy image have a bright central part and therefore 
a larger size. Consequently the faintest surface brightness 
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at the faintest isophotal magnitudes, i.e., Sth — 27.5 mag 
arcsec"^ in Vqqq and 27.0 mag arcsec"'^ in /814, gives the 
surface brightness threshold for object detection. In our 
calculations the above threshold values in the Vgoe and 
1^14 bands are used respectively, although in reality the 
detection was done by the combined Veoe + -^814 image in 
the HDF catalog. 

On the other hand, the surface brightness threshold is 
not clear in the U30Q and B450 bands, because the object 
detection was done without these bands. It should also 
be noted that the isophotal U300 and -B450 magnitudes are 
defined as the flux within the isophotal size in the com- 
bined Vfjoe +-^814 image. Among the objects to which these 
isophotal magnitudes are assigned, those with S/N > 2 
are detected as galaxy images in the U300 and i?45o bands 
(Williams et al. 1996). We see a clear boundary indicated 
by dot-dashed line running from upper-left to lower-right 
in the C/300 and -6450 panels of Fig. ^, and this corre- 
sponds to the line of S/N — 2. It is easy to show that this 
is equivalent to a condition of S + m = const if the noise 
level is proportional to A^/^ (Poisson type noise), where 
A is the isophotal area. We have also confirmed that the 
galaxies in C/300 and -B450 with S/N ^ 2 are actually on 
the dot-dashed line of Fig. || ^. In our calculations the 
galaxies in the C/300 and -8450 bands are detected if they 
are detected in /814 band and furthermore meet the crite- 
rion of S/N > 2 within the isophotal area in the Isu band 
(i.e., those below the dot-dashed line in Fig. |4|). 

2.4. Merger and Number Evolution 

Currently the most popular theory for the structure for- 
mation in the universe is the bottom-up scenario with the 
cold dark matter which dominates the total mass density 
of the universe, in which smaller mass objects form ear- 
lier and then merge into larger objects (e.g., Blumenthal 
et al. 1984). Merging history of dark matter haloes is 
relatively well studied by analytical methods as well as iV- 
body simulations, but merging history of galaxies could 
significantly differ from that of dark matter haloes and is 
poorly known. 

Since a number evolution of galaxies caused by galaxy 
mergers significantly affects the number count of faint 
galaxies, we investigate this effect by using a simple merg- 
ing model in which the luminosity density of galaxies is 
conserved. A common practice for this is to introduce the 
redshift-dependent parameters of Schechter-type luminos- 
ity function of galaxies such as 



(/.*(z) = 0*(O)(l + z)'' 
L*{z) = L*{0){1 + z)-" 



(3) 



We adopt a single value of 77 for all types of galaxies for 
simplicity. 

''This is true for the three WF fields, but not for the 
employed (Williams et al. 1996). The number of galaxies 



In the analysis of this paper, the size of a galaxy is 
crucially important for evaluating the selection effects. 
Merger of galaxies should change their size, and this 
should also be taken into account. The empirical rela- 
tion Te oc Lg may not hold at high redshifts, depend- 
ing on how Te and Lb change during the merger process. 
We assume that the change of Lb and during merger 
processes always satisfy a relation Lb oc r|. (Note that 
this relation is physically different from equation which 
is the relation of galaxies at a fixed time, but describing 
the change of luminosity and size of a test galaxy dur- 
ing merger processes.) If ^ = 2, the surface brightness of 
galaxies is conserved during mergers, and if ^ = 3, the 
luminosity density in each galaxy is conserved. Then the 
luminosity and size of a z = galaxy evolve as 



LB^LBil + z)--^ 
re^re(l + z)-"/« 



(4) 
(5) 



By applying this transformation to the empirical Ts-Lb 
relation at 2: = 0, it is straightforward to give the t^-Lb 
relation as a function of rcdshift: 

re(LB,z) = re(LB,0) X (l + z)-i(f^2-5) , (6) 

which is a generalization of eq. ^. 

The value of ^ depends on the physical process of merg- 
ers. Generally, merger products are expected to become 
more compact than pre-merger progenitors in gas-rich 
mergers because of efficient cooling and dissipation. On 
the other hand, merger products become less compact in 
gas-less mergers because the relative translational energy 
of pre-merger stellar progenitors is converted into the in- 
ternal kinetic energy of a merged stellar system. The for- 
mer corresponds to larger while the latter to smaller f . 
We use ^ = 3 as a standard value and examine the effect 
of changing this value. The assumed conservation of the 
total luminosity density of all galaxies (i.e., (f>*{z)L*{z) = 
constant) will also be discussed when we draw our conclu- 
sion of this paper. 

2.5. Formulations 

In the following we describe the formulations necessary 
to calculate the number count and redshift distribution of 
faint galaxies. We denote A as the observed wavelength, 
and Azj as the wavelength at which the object detection is 
performed. For example, we use A = Ad = Vgoe or /814 for 
counts in the Vgoe and /gi4 bands, whereas A^ — I^u and 
A = C/300 or -6450 for the C/300 and ^450 bands (see §2.3.4). 



The number of galaxies per unit steradian, per unit ap- 
parent magnitude {mx), and per unit redshift is written 

as 



(PN 



dm\dzdQ 



= H{x) 



dzdn ^ 



MMb,z) 



dMB 

dm\ 



(7) 



PC field, because of the different surface brightness threshold 
in the PC field is negligible compared with the WF fields. 
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where x = D\,^ — Umin, -Dad is the isophotal diameter 
of galaxies measured in the band, (PV/dzdVl is the 
comoving volume element which depends on the cosmo- 
logical parameters, ipi is the luminosity function per unit 
Mb with the Schechter parameters given in equation (^, 
and H is the step function [H{x) = 1 and for x > and 
< 0, respectively]. The absolute B magnitude Mb is that 
of the present-day galaxies, which is related with z and 
m\ by i^T-correction and evolutionary [E] correction (see 
equation |lj below). The subscript i denotes the galaxy 
type. The quantity N / dm\dzd^ gives the redshift dis- 
tribution of galaxies, and the integration over z gives the 
number count or the number-magnitude relation. 

We calculate and Mb as a function of to a and z, 
taking into account the selection effects and the photom- 
etry scheme. This is carried out as follows: For a given 
set of Mb and z, by using the assumed luminosity pro- 
file and the Lb-t^ relation, we first calculate the surface 
brightness distribution of a galaxy image in the Ad band. 
Comparing this surface brightness distribution with the 
adopted threshold Sth in the Xjj band, we then calculate 
the isophotal diameter . Given this diameter, we cal- 
culate the isophotal magnitude in the A band. (If the 
aperture magnitude is used in a survey, we should use the 
fixed aperture here. The pseudo-total magnitude is easily 
calculated simply from the total absolute magnitude with- 
out using the surface brightness distribution.) In this way 
we finally obtain to a as a function of AIb and z, or con- 
versely AIb can be related to toa. It is obvious that Z^Ad 
is automatically obtained in this process. In the following 
we give detailed numerical formulations necessary for the 
above calculations. 

Let g{P) be radially symmetric luminosity profile (sur- 
face brightness distribution) of a galaxy, where (3 = r/vg 
is the radius from the center normalized by the effective 
radius of the galaxy. The adopted form of the profile is 
given by 

g{(3) = exp(-a„/3i/") , (8) 

and the integrated profile out to (3 is given by G(/3) = 
27r g{P')l3'df3' . We assume de Vaucouleurs' profile (n = 
4) for elliptical galaxies and an exponential profile (n = 1) 
for spiral galaxies, as mentioned earlier. An effective ra- 
dius Te is defined as the radius within which a half of total 
luminosity is encircled, and by this definition the coeffi- 
cient a„ is given by = 7.67 and ai = 1.68. In order to 
incorporate the effect of observational seeing, we convolve 
this profile function with a Gaussian PSF with dispersion 



f.(0{/0(|)exp 



13^ 



X exp 



2a| 



(9) 



(e.g.. Press et al. 1992). The seeing FWHM in units of 
radian is now related to at as 



seeing FWHM\ ^^(z) 
2^35 / r-e 



(10) 



where dAiz) is the standard angular diameter distance. 

The surface brightness distribution of a galaxy image is 
given by 

S'A(6')[mag arcsec-2] = Ms + (A - B)q + Kx{z) + Ex{z) 
+ 2.5thi(A, z)loge 
+ 51og[re(MB,z)(l + z)Vl0pc] 
+ 26.5721 

-2.51og[.9(/3)/G(oo)] , (11) 

where 9 is the angular radius from the center of a galaxy 
image, K\ is the X-correction, E\ is the i?-correction in- 
cluding internal absorption by dust, and thi is the opti- 
cal depth of intergalactic absorption by HI clouds. The 
present-day color of galaxies of a given type, (A — B)o, 
is calculated by the AB colors based on the present-day 
SEDs of model galaxies which reproduce the observed col- 
ors, and we use -8450 — B = —0.153 to relate the AB mag- 
nitudes to the B magnitudes of the luminosity function. 
The radius parameter (3 is related to 9 as (3 — dA9 jr^, and 
G(oo) = (27i)!7r/a^". The isophotal size can be de- 
rived from solving the equation S\^{D\^ /2) — Sth- Sim- 
ilarly, the apparent magnitude encircled within 9 is given 

by 

mxi9) = Mb + (A - B)o + Kx{\ z) + £;a(A, z) 
+ 2.5rHi(A, z) loge + 5 log[dL/10pc] 
-2.51og[G(/3)/G(oo)] , (12) 

where d^ is the standard luminosity distance, and G{[3) is 
the integrated profile of g{[3). The isophotal, aperture, and 
pseudo-total magnitudes are obtained with 9 = Dy^jl^ 
apcrture/2, and 00, respectively. (Note that G(oo) = 
G(oo).) 

The K and E corrections are calculated from the time- 
dependent models of spectral energy distribution per unit 
wavelength f\{t), after modified to include the dust ab- 
sorption. They can be written as 



i^A(2) = -2.51og 



1 /a/(1+z)(^o) 

(1 + z) fxito) 



Ex{z) = -2. 5 log 



, f f\/{l+z)itz) 



fx/{l + z)ito) 



(13) 
(14) 



where Iq (x) is the modified Bessel function of the first kind 



where is the age of a galaxy at redshift z which was 
formed at Zi;-, and to is its present age. 

With all these prescriptions, we can numerically solve 
Mb for a given set of to a and z, and then the number 
count and redshift distribution of galaxies by means of 
equation (fi\j. 
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3. FAINT GALAXY NUMBER COUNTS 

3.1. Importance of the Observational Selection 

Effects 

Figure ^ shows the predictions of galaxy number count 
in the HDF based on the PLE model of galaxy evolu- 
tion, including the effects of dust and intergalactic absorp- 
tions, and the selection effects. Here we have used a 'stan- 
dard' set of model parameters in this paper: (/i, fig, $7^) = 
(0.7,0.2,0.8), zp = 5, the local luminosity function of 
the SSRS2 survey, and the screen model of dust. [Here, 
h = ifo/(100km/s/Mpc) as usual.] The solid line is the 
total number count of all galaxy types, and the other five 
lines are those of respective galaxy types. In this paper we 
will present many calculations of galaxy counts, changing 
various parameters in order to check the systematic uncer- 
tainties which may affect our conclusions. The summary of 
count calculations is presented in Table || with references 
to the number of figures in this paper. 

In the following we compare these calculated counts with 
the HDF data as well as the ground-based data trans- 
formed into the AB magnitude system. Here we note that 
the HDF bandpass filters for the U300 and Vgoe are signifi- 
cantly different from those for the U, V, and R used in the 
ground-based observations. We have corrected the magni- 
tudes of such ground-based data, making use of the central 
wavelength of bandpass filters, the present-day SEDs of 
different galaxy types, and the relative proportions among 
the galaxy types. The HDF counts are those of isophotal 
magnitudes, while published total-magnitude counts are 
used for other ground-based data for which the selection 
effects are not as important as for the HDF data. 

Figure || shows the effect of cosmological parameters 
and the selection effects in the predictions of galaxy num- 
ber count. The predictions with the selection effects are 
presented by solid lines, while those without the selec- 
tion effects by dashed lines. In either cases, the three 
lines from bottom to top correspond to the EdS uni- 
verse with (h, r^o, JIa) ~ (0-5, 1, 0), an open universe with 
{h, fio) ^^a) = (0.6, 0.2, 0), and a A-dominated flat universe 
with {h,no,^A) = (0.7,0.2,0.8). 

The dotted line is a prediction without the selection ef- 
fects in an open universe with CIq — 0.1. The prescrip- 
tions used here are the same as those in the prediction 
by Pozzetti et al. (1998) which they found to agree with 
the observed HDF counts. Our prediction shown by this 
dotted line also gives a good fit to the HDF data, and we 
have confirmed Pozzetti et al.'s result where no selection 
effects are taken into account. This suggests that the dif- 
ference between the PLE models of different authors is not 
significant in the number count predictions. 

However, this figure clearly demonstrates the impor- 
tance of the selection effects in comparison between the 
predicted and observed number counts of galaxies. The 



difference between the predicted counts with and without 
the selection effects attains up to a factor of about 4 at the 
faintest magnitudes. As a result, the observed HDF counts 
in excess of the PLE predictions with the selection effects 
are much larger than previously considered. It is striking 
that such PLE predictions seem to be short of the observed 
counts in all the passbands, even in the A-dominated flat 
universe with JIa — 0.8 in which the number count is close 
to the maximum. In the EdS universe, such a predicted 
deficit in the PLE model attains up to a factor of more 
than 10. This large excessive number of HDF galaxies may 
suggest the number evolution of galaxies, but it depends 
heavily on the cosmology how much the number evolution 
is required to explain the HDF data. 

3.2. Dependence on Model Parameters 

Before we investigate the effect of number evolution 
quantitatively, it is necessary to see the uncertainties in 
the PLE predictions. Figure ^ shows how the PLE pre- 
dictions depend on models of luminosity evolution and ab- 
sorptions. The solid line is the standard model presented 
in Fig. H The dotted line is the same as the standard 
model but with no luminosity evolution. The dot-dashed 
line is the prediction where the updated luminosity evolu- 
tion model of Kobayashi et al. (1999) is used rather than 
the standard model of Arimoto & Yoshii (1987) and Ari- 
moto, Yoshii, & Takahara (1992). The short dashed line 
is the prediction where the slab distribution of dust is as- 
sumed rather than the screen by dust. The long-dashed 
line is the prediction where no intergalactic absorption is 
taken into account. This line almost overlaps with the 
solid line and is not visible except in the U300 band. 

Figure || shows the effect of changing the galaxy forma- 
tion epoch zp and local luminosity function. The solid 
line is the standard model presented in Fig. ||. The short- 
and long-dot-dashed lines are the predictions with zp — 3 
and 10, respectively, rather than zp = 5 in the standard 
model. The short- and long-dashed lines are the predic- 
tions where the luminosity function of the Stromlo-APM 
and the CfA redshift surveys are used, respectively, with 
zp = 5. 

Figure ^ shows how the uncertainty in the size- 
luminosity relation affects the predictions of galaxy num- 
ber count. The solid line is the standard model presented 
in Fig. 1^, while the dashed line is the prediction when 
the selection effects are ignored, as shown in Fig. ||. The 
short- and long-dot-dashed lines are the counts when the 
r^-Ls relation is shifted by ± 1 cr in A (log re). The dotted 
line is the prediction using the GCE sequence instead of 
the standard GDE sequence as the size-luminosity relation 
of elhptical galaxies. (See § p.3.1 for detail.) 

From these three figures we understand a range of the 
systematic uncertainties in the PLE predictions, which 
turns out to be not large enough to save the EdS uni- 
verse nor an open universe with fio > 0.2. These uncer- 



8 



tainties are significant at the faintest magnitudes where 
the N-m relation starts to turn over. However, we point 
out that the effect of cosmological parameters becomes 
apparent at brighter magnitudes where the uncertainties 
remain much less significant. In fact, we will show that the 
slope of the N-m relation at i?45o — 22-26 can be used to 
discriminate between the effects of cosmological parame- 
ters and galaxy number evolution. In order to demon- 
strate the above statement quantitatively, the systematic 
model uncertainties at m= 25 and 28 are summarized in 
Table ^ The uncertainty at m = 25, to be compared with 
the effect of cosmological parameters, is dominated by the 
dust distribution model, but we note that the estimated 
change from a standard screen model to the slab model is 
somewhat overestimated because the slab model is clearly 
inconsistent as a model of dust distribution in starburst 



galaxies, as mentioned in {2.2 



3.3. Mergers and Number Evolution 

The PLE predictions fall considerably short of the ob- 
served HDF counts in the EdS universe and an open uni- 
verse, and it is still the case even in a A-dominated fiat 
universe. Here, by using a simple model of mergers intro- 



duced in §2.4, we investigate whether the number evolution 
explains the large number of faint HDF galaxies. Figure 



10 shows the effect of introducing such a number evolution 
model in a A-dominated fiat universe. The solid line is the 
prediction with the merger parameters of (ry,^) = (1,3), 
while the dotted line is the standard PLE prediction with 
no number evolution (Fig. |^). The short- and long-dot- 
dashed lines are the predictions with (77,^) — (1,2) and 
(1, 4), respectively, showing that the effect of changing 
^ is not significant. This result indicates that a modest 
number evolution with ~ 1 [</)* cx (1 -I- z)^\ is sufficient 
to explain the observed HDF counts in a A-dominated fiat 
universe with JIa = 0.8, and an even stronger number 
evolution with 77 > 1 is rejected by the data. 

Next we consider the EdS universe where the PLE count 
prediction is by more than one order of magnitudes smaller 
than those observed in the HDF (Fig. |ll|). The predic- 
tions with = 2, 3, 4, and 5 are shown by four solid lines in 
order from bottom to top, with a fixed value of ^ = 3. The 
dotted line is the PLE prediction without number evolu- 
tion. This result indicates that a strong number evolution 
with > 3-4 is necessary to explain the HDF counts in 
the EdS universe. We note that, while the strong number 
evolution explains the counts at the faintest magnitudes, 
it fails to explain the overall shape or slope of the N-m re- 
lation. This failure is clearly seen in the -6450 band, where 
the strong number evolution makes the count slope less 
steep and deviate from the observations most prominently 
at ^450 = 22-26. This argument is quite robust, because 
there should be an upper bound in the steepness of the 
N-m slope (dlogN/dm < 0.4), when the total luminos- 
ity density of all galaxies is conserved during the merger 



process. In a more realistic case such as gas-rich merg- 
ers inducing starbursts, an even flatter slope is predicted, 
because galaxies before mergers are always fainter than 
those in the case of luminosity-density conservation. If 
one tries to explain the observed steep slope by mergers, 
it is necessary to contribe a merger process where galaxies 
before mergers are always brighter, in other words, the lu- 
minosity density of galaxies increases if the merger process 
is traced backwards — which we consider quite unrealistic. 
Therefore, we conclude that a strong number evolution in 
the EdS universe is unlikely to explain the observed counts 
over the whole range of apparent magnitudes. 

Figure ^ shows the effect of introducing a number evo- 
lution model in an open universe with Qq = 0.2. The lines 
in this figure have the same meanings as in Fig. |ll], but 
the four solid lines are the predictions with 77 = 1, 2, 3, 
and 4 in order from bottom to top. In this open universe, 
a number evolution with > 2 is necessary to explain the 
faintest counts, but again it can not explain the steep slope 
of the N-m relation in the and Vgoe bands. This indi- 
cates that an open universe is also difficult to explain the 
HDF counts if > 0.2. However, a lower-density open 
universe, for example, with ~ 0.1 might give a similar 
result with a A-dominated fiat universe with Oq = 0.2 (see 
Fig. ^), and such an open universe could also explain the 
HDF counts if a modest number evolution of galaxies is 
taken into account. 



4. PHOTOMETRIC REDSHIFT DISTRIBUTION 

Although not as reliable as spectroscopic redshifts, pho- 
tometric redshifts of galaxies are useful for statistical stud- 
ies of high-redshift galaxies. Several groups have pub- 
lished catalogs of photometric redshifts for the HDF galax- 
ies (e.g., Sawicki, Lin, & Yee 1997; Wang, Bahcall, & 
Turner 1998; Fernandez-Soto, Lanzetta, & Yahil 1999). 
Here we compare our theoretical model with the photo- 
metric redshift distributions reaching /814 — 28 derived by 
Fernandez-Soto et al. (1999), which utilizes not only the 
optical photometry of the HDF but also the information 
of the near-infrared J, H, and K bands. We have also 
compared our model with the other two catalogs of photo- 
metric redshifts by Sawicki et al. (1997) and Wang et al. 
(1998), and confirmed that the following result is hardly 
changed. 

Figure |l^ shows the observed redshift distribution in 
three /814 magnitude ranges. The model curves are 
our predicted redshift distributions of galaxies in a A- 
dominated fiat universe with Q\ = 0.8. The area over 
which the models are calculated is chosen to coincide 
with the sky area covered by the observational analysis of 
Fernandez-Soto et al. (1999): 5.31 arcmin^ for Isu < 26 
and 3.92 arcmin^ for /gi4 > 26. Therefore, not only the 
shape but also the normalization of the predicted redshift 
distributions can be compared directly with the data. The 
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solid and dashed lines are the models with = (0,3) 

and (1,3), respectively, with the selection effects taken into 
account. The dot-dashed line is the same as the dashed 
line with (r;,^) = (1,3), but no selection effects are taken 
into account. It is clear that the selection effects give a 
bias against high-redshift galaxies, and this selection bias 
is significant especially at the faintest magnitudes. It is in- 
evitable to include these effects when one uses the redshift 
distribution as a probe of number evolution of galaxies. 
Comparison with the data shows that a modest merger 
model with 77 1 in a A-dominated flat universe gives a 
reasonable fit to the photometric redshift distribution as 
well as galaxy counts, provided that the selection effects 
are properly taken into account. 

Figure |lj compares the observed redshift distribution 
with the predictions in the EdS universe. The solid line 
is the prediction without number evolution, while other 
lines are those with 77 = 3, 4, and 5 with a fixed value of 
^=3. The selection effects are taken into account in all 
curves. A strong number evolution predicts that most of 
galaxies have lower redshifts of 2: < 1, deviating signifi- 
cantly from the observed distribution for 23 < /814 < 26. 
If the assumption of conserved luminosity density is re- 
laxed, the predicted distribution becomes peaked at an 
even lower redshift, because pre-merger galaxies in more 
realistic merger models are fainter than expected from the 
conserved luminosity density, as discussed in the previous 
section. This gives another argument that the EdS uni- 
verse can not explain the observed number of HDF galax- 
ies even if a strong merger is invoked, in addition to the 
argument in the previous section against the EdS universe 
based on the slope of the observed N-m relation. 

Figure |l^ is similar to Fig. but for an open uni- 
verse with rio = 0.2. The predicted redshift distribu- 
tion is still peaked at a lower redshift compared with the 
observed distribution. This discrepancy is significant for 
23 < /814 < 26, but not as serious as in the EdS universe. 

5. DISCUSSION 

In this paper we have shown that a strong number evolu- 
tion with rj > 3-4 is necessary to explain the HDF counts 
at the faintest magnitudes in the EdS universe. The photo- 
metric redshift distribution of HDF galaxies suggests that 
a significant number of the faintest galaxies are at z < 1, 
and hence there must be a strong number evolution al- 
ready operated at z < 1 to increase the number of galaxies 
by factor of about 10. Therefore, further argument for or 
against the EdS universe with strong number evolution of 
galaxies can be made from the observational constraints 
based on spectroscopic redshift surveys at z < 1. 

Totani & Yoshii (1998) argued that such a strong num- 
ber evolution at z < 1 is clearly inconsistent with the 
spectroscopic catalogs of galaxies reaching z ~ 1, at least 
for giant galaxies with L ~ L*. This is based on a F/Vmax 
test for the galaxies in the Canada-France Redshift Sur- 



vey (CFRS, Lilly et al. 1995), and an important finding 
is that the PLE model is not inconsistent with the spec- 
troscopic data, giving a constraint on number evolution as 
77 = 1.8 ± 0.7, 1.1 ± 0.7, and 0.5 ± 0.7 for the EdS universe, 
an open l^o = 0.2 universe, and a flat J^a = 0.8 universe, 
respectively. Although the F/14nax test favors a larger 77 
in the EdS universe, this value seems smaller than that 
required to explain the HDF counts at the faintest mag- 
nitudes. These results have been confirmed for elliptical 
galaxies by Shade et al. (1999), in which they found that 
the population of massive early-type galaxies was largely 
in place by z 1. 

An independent constraint on galaxy mergers at z < 1 
comes from statistical studies of merging galaxies inferred 
from high-resolution images. Recently, Le Fevre et al. 
(2000) derived the evolution of merger rate from the HST 
images of the CFRS and LDSS galaxies. Their result 
suggests that L* galaxies on the average have undergone 
about one merger event from z = 1 to 0, which corresponds 
to 77 1 and hence this is consistent with the number evo- 
lution in a A-dominated flat universe suggested by this 
paper. 

Pozzetti et al. (1998) claimed that the PLE model can 
not explain all the observed data, although it well explains 
the HDF galaxy counts. The major discrepancies between 
the PLE model and the observed data were found in the 
evolution of the luminosity density in the universe. That 
is, the observed luminosity density increases more steeply 
to z ^ 1 than the PLE model prediction, and on the other 
hand, the PLE model predicts too high UV luminosity 
density at 3.5 < z < 4.5 compared with the observation, 
because of intense starbursts in elliptical galaxies. How- 
ever, Totani, Yoshii, & Sato (1997) had already pointed 
out that the observed steep evolution to z ~ 1 is explained 
in a A-dominated flat universe. It has already been argued 
that the PLE model is consistent with the observed lumi- 
nosity density evolution if initial starbursts in high-redshift 
elliptical galaxies are obscured or not existent. In fact, our 
calculation of the redshift distribution in this paper shows 
that initial starbursts are not detected at z > 3, because 
of the dust obscuration and the selection effects. There- 
fore, the major problems of the PLE model claimed by 
Pozzetti et al. (1998) are resolved, and the PLE model in 
a A-dominated flat universe gives a reasonable fit to the 
observed HDF data, allowing only for a modest number 
evolution of galaxies with 77 < 1. 

6. SUMMARY AND CONCLUSIONS 

We have modeled the number count and redshift dis- 
tribution of faint HDF galaxies, with the observational 
selection effects properly taken into account. As a conse- 
quence of the selection effects in the theoretical modeling, 
predicted counts from the PLE model are smaller than 
previously considered, and they are more than 10 times 
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smaller than the observed HDF counts at the faintest mag- 
nitudes in the EdS universe. A strong number evolution 
with 77 > 3-4 under the assumption of conserved luminos- 
ity density is required to explain the faintest counts in this 
EdS universe, when the number evolution is parametrized 
as ^* oc {1 + z)^ and L* oc (1 -I- z)~^. However, such a 
strong number evolution is not consistent with the overall 
N-m slope or the photometric redshift distribution. These 
discrepancies become even worse when one considers a 
more realistic merger process, i.e., enhanced star forma- 
tion following by gas-rich mergers. In addition, such a 
strong evolution is rejected at least for average L* galax- 
ies at z < 1 from the data of spectroscopic redshift surveys 
(Totani & Yoshii 1998; Shade et al. 1999; Le Fevre et al. 
2000). Therefore, we conclude that it is almost impossible 
to explain the HDF galaxies in the EdS universe, unless 
we invoke ultra-exotic galaxy populations such as galaxies 
forming only massive stars at high redshifts to escape from 
local galaxy surveys due to the complete lack of long-lived 
stars. 

The present work revitalizes the practice of using faint 
number counts as an important cosmological test, which 
gives one of the arguments against the EdS universe by 
its outstanding statistics compared with other cosmologi- 
cal tests. An open universe with Qq > 0.2 does not fit to 
the HDF data either, for the similar reasons for rejecting 
the EdS universe. An open universe with ~ 0.1 might 
be consistent with the HDF data, but such a low value of 
Qq ~ 0.1 would not be reconciled with other constraints 
on r^o , such as the baryon-gas to dark-matter mass ratio in 
poor clusters of galaxies combined with the standard big- 
bang nucleosynthesis prediction of baryon mass density 
in the universe (e.g., Pedersen, Yoshii, & Sommer-Larsen 
1997). 

We have extensively checked systematic uncertainties in 
our theoretical modeling of galaxy formation and evolu- 
tion, and found that they are unlikely to resolve the above 
discrepancies emerged in the EdS universe and also in an 
open universe. On the other hand, such discrepancies are 
naturally resolved if we invoke a A-dominated flat uni- 
verse. This suggests that the existence of the cosmological 
constant or an exotic form of the vacuum energy density 
of the universe which is now accelerating the expansion of 
the universe. 

The PLE model in a A-dominated flat universe with 
fio ^ 0.2 gives a reasonable fit to the HDF data, and 
a modest number evolution with ry < 1 is also suggested 
by the HDF counts at the faintest magnitudes. It should 



be noted that this number evolution does not necessar- 
ily mean mergers of galaxies, but may suggest strongly 
clumpy star-forming regions within an individual galaxy 
system becoming visible at high redshifts (CoUey et al. 
1996; Bunker, Spinrad, & Thompson 1999). 

On the other hand, it is interesting to note that this 
indication of mild number evolution is consistent with the 
merger rate evolution of L* galaxies at 2; < 1 recently in- 
ferred from a high-resolution image study for galaxies in 
the CFRS survey (Le Fevre et al. 2000). This result is 
consistent with some models of galaxy formation based on 
the hierarchical structure formation in the CDM universe 
(Le Fevre ct al. 2000), although there arc considerable 
uncertainties in the theoretical calculations for the merg- 
ing history of baryonic component. A stronger number 
evolution with t] > 1 is, however, strongly disfavored by 
the observed HDF galaxy counts. This will give an im- 
portant constraint when galaxy formation is modeled in 
the framework of the structure formation in a cold dark 
matter universe. 

Inclusion of the selection effects in this paper leads 
to a considerably different result from previous studies 
on galaxy number count and redshift distribution. This 
means that any cosmological interpretations will be seri- 
ously misled if the selection effects are ignored. All fu- 
ture studies related to the detection and statistics of high- 
redshift galaxies should take into account these effects. 
The selection effects give a bias against high-redshift galax- 
ies, reducing a problem of overprediction of such galaxies 
by the PLE model, which has been claimed by several stud- 
ies ignoring the selection effects (e.g., Ellis 1997). In fact, 
we have shown that the PLE model is in overall agree- 
ment with the HDF galaxies, even if a modest number 
evolution of galaxies (ry < 1) may be required. A strong 
number evolution, however, predicts too small a number 
of high-rodshift galaxies to be consistent with the photo- 
metric redshift distribution of the HDF galaxies. 
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Table 1 

schechter parameters adopted for local luminosity functions 



Survey 


Type 


(j)* [Mpc-^] 


a 


M* 


References 


SSRS2 


E/SO 


4.4 X 10"^ 


-1.00 


-19.37 


1 




Sab-Scd 


8.0 X 10^3 


-1.11 


-19.43 






Sdm 


0.2 X 10-3 


-1.81 


-19.78 




Stromlo-APM 


All 


1.4 X 10-2 


-0.97 


-19.50 


2 


CfA 


E/SO 


4.6 X 10-3 


-0.88 


-19.23 


3 




Sab 


3.8 X 10-3 


-1.00 


-18.99 






She 


6.6 X 10-3 


-0.44 


-19.06 






Scd 


2.0 X 10-3 


-1.00 


-19.34 






Sdm 


0.8 X 10-3 


-1.20 


-19.29 





NOTE. — Hq = 100 km/s/Mpc. SSRS2 is adopted as a standard in our analysis. The value of Mg is given for the B 
band, which is related to i?45o as i?45o = B — 0.153. 

REFERENCES.— (1) Marzke et al. (1998), (2) Loveday et al. (1992), (3) Huchra et al. (1983). See also Efstathiou, Ellis, 
& Peterson (1988) and Yoshii & Takahara (1989). 



Table 2 

Figure Identification for Galaxy Count Calculations 



Examined Effects 



Figure No. 



Galaxy morphological types 

Selection effects 

Cosmological parameters 

Galaxy luminosity evolution 

Absorption (dust and intergalactic HI clouds) 

Local luminosity function 

Formation redshift zp 

Dispersion in the Ls-Ve relation 

Merging of galaxies 



3 



io|^y 



Table 3 

Sensitivity of the Predicted N-m Relation to the Change of Input Parameters 



Veoe 



# 


Change 


A log Af25mag 


A log Af28mag 


A log Af25mag 


A log Af28mag 


A log iV25mag 


AlogAr28 


1 


Cosmology: A — > EdS/open" 


-0.227/-0.126 


-0.513/-0.232 


-0.446/-0.235 


-0.610/-0.278 


-0.493/-0.254 


-0.601/-0 


2 


Evolution: AY, AYT^KTN*- 


+0.048 


-0.026 


-0.027 


-0.086 


-0.056 


-0.017 


3 


Dust absorption: screen— >slab 


-0.085 


-0.166 


+0.062 


+0.101 


+0.137 


+0.13( 


4 


HI absorption: on^ofF 


+0.094 


+0.138 


+0.020 


+0.009 


+0.005 


+0.00^ 


5 


Local LF: SSRS2^Stromlo-APM/CfA'= 


-0.119/+0.066 


-0.283/-0.225 


-0.071/+0.041 


-0.218/-0.128 


-0.042/+0.022 


-0.187/-0 


6 


Formation epoch zp: 5 ^ 3/10 


+0.054/-0.040 


-0.101/+0.036 


-0.070/-0.062 


-0.151/+0.075 


-0.119/-0.025 


-0.151/+0 


7 


re-Lp relation: +la/—la in A(logre)'^ 


-0.088/+0.045 


-0.159/+0.103 


-0.150/+0.069 


-0.160/+0.098 


-0.104/+0.050 


-0.147/+0 




Total systematic uncertainty*^ 


±0.157 


±0.339 


±0.155 


±0.277 


±0.185 


±0.27! 



"A: (r2o,!^A) = (0.2, 0.8), EdS: (1, 0), open: (0.2, 0) 

''AY: Arimoto & Yoshii (1987), AYT: Arimoto, Yoshii, & Takahara (1992), KTN: Kobayashi, Tsujimoto, & Nomoto (1999). 
'^See Table F] for detail. 
''See for detail. 

■^Quadratic sum of the rows 2, 3, 5-7. A mean value of the two numbers in the rows 5—7 is used in the sum. 



NOTE. — The prescriptions of a standard model in our analysis include a A cosmology, AY-AYT evolution model, the local LF of SSRS2, 

Zp = 5, the screen model of dust absorption, a nd the intergalactic HI absorption. 

[Complete version of this table is available at http://th.nao. ac.jp~totani/imagcs/paper/ty2000-table3.ps ] 
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103 lo" 10= 10* IQi* 

Wavelength [A] 



Fig . 1 . — Spectral energy distributions (SEDs) of various morphological types of galaxies at epochs of 1 , 5 , and 1 5 Gyrs 

(from the top to the bottom) after their formation. For the purpose of clarity, the SEDs at 1 and 5 Gyrs are artificially 
multiplied by a factor of 10'* and 10^, respectively. The screen model of dust extinction is used as our standard and 
is incorporated in the calculations shown by the solid line. The dashed and dotted lines are those from the slab dust 
model and the no-extinction model, respectively. The SEDs at 15 Gyrs are made to agree with the observed SEDs in 
the present-day galaxies. 




2000 4000 6000 8000 10000 12000 14000 
Wavelength [A] 



Fig. 2. — Optical depth of intcrgalactic absorption by intervening HI clouds as a function of observed wavelength for 
various source redshifts of z = 1, 2, 3,..., 10 (from left to right). Shown are the calculations by Yoshii & Peterson (1994), 
with the Doppler 6-parameter of 6 = 20 km/s in the HI clouds (Spitzer 1978). 
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Fig. 3. — Size-luminosity relations for local elliptical and spiral galaxies. The data are from Bender et al. (1992) 
for elliptical galaxies and Impey et al. (1996) for spiral galaxies. For elliptical galaxies, the filled circles are for giant 
ellipticals, the open circles for intermediate ellipticals, the open squares for bright dwarf ellipticals, the filled triangles 
for compact ellipticals, and the open pentagons for dwarf spheroidals. The solid lines are the best-fit relations for the 
giant-dwarf elliptical sequence (GDE) and spiral galaxies, while the dot-dashed line is the best-fit to the giant-compact 
elliptical sequence (GCE). The dashed lines are the relation shifted by standard deviation in A(logre) from the solid 
lines. 
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Fig. 4. — Isophotal AB magnitude versus average surface brightness of the HDF galaxies (Williams et al. 1996). 
The dashed line is the limit of minimum area of a galaxy image A = 0.04 arcsec"^. The dot-dashed lines in the Veoa 

and /814 bands arc the adopted isophotal thresholds {Sth = 27.5 and 27.0 mag arcscc"^ for the Veoa and /8i4 bands, 
respectively) . The dot-dashed lines in the t/300 and bands correspond to the detection threshold by signal to noise 
ratio of S/N > 2 (see text). 
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Fig. 5. — Faint galaxy number counts predicted in the four bands of the HDF. The selection effects under the HDF 
observation condition are included. The model presented here is our 'standard' PLE model with the standard set of 
model parameters: {h, Clo, ^a) = (0.7, 0.2, 0.8), zf = 5, the local luminosity function of the SSRS2 survey, and the screen 
model of dust extinction. Intergalactic absorption by HI clouds is also taken into account. The solid line is the total 
counts of all five types of galaxies, while the other lines are for individual types of E/SO (dotted). Sab (short-dashed), 
Sbc (long-dashed), Scd (dot-short-dashed), and Sdm (dot-long-dashed). 
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Selection Effects in Different Cosmological Models 




20 25 30 20 25 30 
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Fig. 6. — Comparison of the PLE model with the observed HDF counts as weU as ground-based brighter counts. The 
top solid line is the standard PLE model shown in Fig. ||in a A-dominated flat universe with {h, ^Iq, JIa) = (0.7, 0.2, 0.8), 
while the middle and bottom lines are for the different cosmologies of an open universe with (0.6, 0.2, 0.0) and an EdS 
universe with (0.5, 1, 0), respectively. The three dashed lines from top to bottom are the same as the solid lines, except 
that the observational selection effects are not taken into account. The dotted line is a prediction in an open universe 
with (/i, r^Oj^^A) = (0.6,0.1,0.0) and no selection effects, which are the same prescription with the previous study by 
Pozzetti et al. (1998). 
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Fig. 7. — Comparison of the PLE model with the observed HDF counts as well as ground-based brighter counts. 
The solid line is the standard PLE model shown in Fig. ||. The other lines are the same but for the case of no 
luminosity evolution (dotted), slab model dust (short-dashed), no intergalactic absorption by HI clouds (long-dashed). 
The dot-dashed hne is for the PLE model with the updated luminosity evolution model of Kobayashi et al. (1999). 
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Effects of Different Local LF and Formation 
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Fig. 8. — Comparison of the PLE model with the observed HDF counts as well as ground-based brighter counts. The 
solid line is the standard PLE model shown in Fig. |^. The other lines are the same but for the models with zp — 3 
(short-dot-dashed) and 10 (long-dot-dashed), and with different local luminosity functions of the Stromlo-APM survey 
(short-dashed) and the CfA survey (long-dashed). 
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Fig. 9. — Comparison of the PLE model with the observed HDF counts as well as ground-based brighter counts. The 
solid line is the standard PLE model shown in Fig. |^. The other lines are the same but for models without the selection 
effects (dashed) and with the shifted size-luminosity relation by +la (short-dot-dashed) and —la (long-dot-dashed) 
deviation in A (log re). The dotted line is the PLE model where the GCE sequence of elliptical galaxies is used instead 
of the standard GDE sequence. See §2.3.1 for detail. 
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Fig. 10. — The effect of galaxy number evolution on galaxy counts based on the luminosity-density conserving mergers 
in a A-dominated universe with {h, ilo, Ha) — (0.7, 0.2, 0.8). The dotted line is the standard PLE model without number 
evolution in this universe shown in Fig. |^ as well as Fig. ^ The short-dot-dashed, solid, and long-dot-dashed lines are 
the models with the merger parameters of (??,^) =(1, 2), (1, 3), and (1, 4), respectively. 
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Fig. 11. — The effect of galaxy number evolution on galaxy counts based on the luminosity-density conserving mergers 
in the Einstein-de Sitter (EdS) universe with {h, Hq, Ha) — (0.5, 1, 0). The dotted line is the standard PLE model without 
number evolution in the EdS universe shown in Fig. |^. The four solid lines are the models with the merger parameters 
of r] = 2, 3, 4, and 5 (from the bottom to the top), where ^ is fixed to 3. 
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Merging in an Open Universe 
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Fig. 12. — The effect of galaxy number evolution on galaxy counts based on luminosity-density conserving mergers in 
an open universe with (/i, r2o,i^A) — (0.6,0.2,0). The dotted line is the standard PLE model without number evolution 
in an open universe shown in Fig. ^ The four solid lines are the models with the merger parameters of rj = 1, 2, 3, and 
4 (from the bottom to the top), where ^ is fixed to 3. 
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Photo-z Distribution in a A Universe 




Redshift z 

Fig. 13. — The comparison with the photometric redshift distribution derived by Fernandez-Soto, Lanzetta, & Yahil 
(1999) in a A-dominated flat universe with {h, fto, J7a) = (0.7, 0.2, 0.8). The sohd hne is the standard PLE model without 
number evolution, with the same model parameters as in Fig. |^. The dashed line is the merger model with the merger 
parameters of (77, C) = (1;3), which reproduces the observed HDF number counts as well. The dot-dashed line is the 
same with the dashed line, but the selection effects are not taken into account. 



Photo-z Distribution in the EdS Universe 
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Fig. 14. — The comparison with the photometric redshift distributions derived by Fernandez-Soto, Lanzetta, & Yahil 
(1999) in the Einstein-de Sitter universe with {h,no,n\) — (0.5, 1,0). The solid line is the PLE model without number 
evolution. The short-dashed, long-dashed, and dot-dashed lines are the merger models with r/ = 3, 4, and 5, respectively, 
where ^ is fixed to 3. The selection effects are taken into account in all these models. 
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Photo-z Distribution in an Open Universe 




Redshift z 

Fig. 15. — The comparison with the photometric redshift distributions derived by Fernandez-Soto, Lanzetta, & Yahil 
(1999). The same as Fig. ^ but for an open universe with (ft,, JIqi^^a) — (0.6,0.2,0). The sohd hne is the PLE model 
without number evolution. The short-dashed, long-dashed, and dot-dashed lines arc the merger models with r; = 1, 2, 
and 3, respectively, where f is fixed to 3. The selection effects are taken into account in all these models. 



